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Cooperative effects of Coulomb and electron-phonon interactions in the 
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We study the electronic states and the superconductivity in the two-dimensional 16-band d-p 
model coupled with Aig, Big and Eg local phonons and obtain the rich phase diagram including the 
magnetic, charge and orbital ordered phases on the parameter plane of the Coulomb and electron- 
phonon interactions. When the electron-phonon interaction is dominant, the charge fluctuations 
induce the s+_|.-wave superconductivity, while when the Coulomb interaction is dominant, the mag- 
netic fluctuations induce the s±-wave superconductivity. Remarkably, the orbital fluctuations are 
enhanced due to the cooperative effects of the Coulomb and electron-phonon interactions and induce 
the s-|-+-wave and the nodal s±-wave superconductivities. 



I. INTRODUCTION 

The recently discovered iron-based superconductors^'^ 
RFePnOi-a^Fj: (R=Rare Earth, Pn=As, P) with a 
transition temperature Tc exceeding 50K— "— have at- 
tracted much attention. At present, there are fol- 
lowing four families of the iron-based superconductors: 
RFeAsO with ZrCuSiAs-type structure (1111 system), 
BaFe2As2 with ThCr2Si2-type structure (122 system)^i2,, 
LiFeAs and NaFeAs with PbFCl-type structure (111 
system)iS-iS and Fe(Se,Te) (11 system)^.^^. These sys- 
tems have similar conducting Fe-pnictogen (-chalcogen) 
planes and the resulting electronic structures predicted 
by the first principle calculations are similar to those 
for each familiesi^""— . The energy bands near the Fermi 
level mainly constructed by the Fe-3d orbitals are heav- 
ily entangled and there are two or three concentric hole 
Fermi surfaces (FSs) around the F-point [k = (0, 0)] 
and the elliptical electron FSs around the A/-point [fc = 
(tt, tt)]. These features are observed by the angle re- 
solved photoemission spectroscopy (ARPES) in several 
compounds?-:^"^ . Despite the similarities of the electronic 
structures for the four families, it seems that the details of 
the gap structures are different from system to system as 
mentioned below and the pairing state together with the 
mechanism of the superconductivity for the iron-based 
superconductors is still controversial. 

As for the 1111 system, the F nondoped compound 
LaFeAsO exhibits the structural transition from tetrag- 
onal (P4/nmm) to orthorhombic (Cmma) phase at a 
transition temperature T =155K and stripe-type anti- 
ferromagnetic (AFM) order at T = 134K with a mag- 
netic moment ^ 0.36/i5^^ at low temperature. With in- 
creasing F doping, the system becomes metallic and the 
AFM order disappears^, and then, the superconductivity 
emerges for x ^ 0.1 with Tc ~ 26K. Rare-earth substitu- 
tion compounds exhibit superconducting transition with 
higher T^—. The NMR Knight shift measurements re- 
vealed that the superconductivity of the systems is the 
spin-singlet pairing^'*-^^. Fully-gapped superconducting 
states have been predicted by various experiments such 
as the penetration deptb^S, the specific heat^, and the 



impurity effect on T^-^. In contrast to the above men- 
tioned experiments, the NMR relaxation rate shows the 
power low behavior l/Ti ex T^ below T^M-, suggesting 
the nodal or highly anisotropic gap structure. The other 
NMR measurements^i, however, revealed 1/Ti ex T^ be- 
low Tc and there is still controversy. 

The parent compound of the 122 system BaFe2As2 
shows the structural transition from tetragonal 
(M/mmm) to orthorhombic (Fmnim) phase and 
the stripe-type AFM order simultaneously at a tran- 
sition temperature T ^lAQY^^^, where the magnetic 
moment is about 0.87piB at low temperature^. Both 
the electron and hole doping by the substitution Co for 
Fe and K for Ba induce the superconductivitySi^. The 
T^-dependence of 1/Ti— , the exponential behavior of 
the penetration depth^S. and the ARPES^ii^ada suggest 
the fully-gapped superconductivity. The chemical 
pressure by substituting P for As in BaFe2(Asi_;rPa;)2 
also leads to the superconductivity with Tc up to 30K^«^ 
where the specific heat, the penetration depth, the 
thermal conductivity and NMR 1/Ti imply the nodal or 
highly anisotropic gap structures4^i^. 

In 11 system, FeTe shows the another type of the AFM 
order with ordering vector q — (tt, 0) ^4,45 ^ where the 
magnetic moment is about 2.03/ifl'*— i^. On the other 
hand, FeSe does not exhibit the magnetic order but su- 
perconducting transition at T ~ 8Ki^. The thermal 
conductivity^^ suggests the fully-gapped superconduct- 
ing state. 

Theoretically, Mazin et al. suggested that the fully- 
gapped s-wave pairing whose order parameter changes 
its sign between the hole FSs and the electron FSs {s±- 
wave pairing) is favored due to the stripe-type AFM 
spin fluctuations^. According to the weak coupling ap- 
proaches based on multi-orbital Hubbard models^"—, 
and those based on the d-p model§2r— , the s±-wave pair- 
ing seems to be the promising candidate for the pairing 
state in the iron-based superconductors. It is shown that 
the s±-wave pairing is realized also in the strong cou- 
pling region by the mean field study based on the f- Ji- J2 
model^ and the exact diagonalization study based on the 
one-dimensional two-band Hubbard model^^. The s±- 



wave state mediated by the spin fluctuations seems to be 
consistent with many experiments. However, the theoret- 
ical analysis of the nonmagnetic impurity eflfects based on 
the 5-band Hubbard model shows that the s±-wave state 
is very fragile against nonmagnetic impurities^. This 
is in contradiction to the experimental results that the 
superconductivity for the iron-based superconductors is 
robust against nonmagnetic impurities. Therefore, the 
fully-gapped s-wave state without sign reversing (s++- 
wave state) is considered to be another promising candi- 
date for the pairing state in the iron-based superconduc- 
tors. 

In the previous papers^"— i^, we have investigated 
the electronic states of the Fe2As2 plane in iron-based 
superconductors on the basis of the two-dimensional 16- 
band d-p model which includes the Coulomb interaction 
on a Fe site: the intra- and inter-orbital direct terms U 
and U' , the Hund's coupling J and the pair-transfer J'. 
Using the random phase approximation (RPA), we have 
found that, for U > U', the s±-wave superconductivity 
is realized due to the spin fluctuations with q ^ (7r,7r), 
while for U < U' , the s++-wave state is realized due 
to the orbital fluctuations with q = (0, 0)^^. In addi- 
tion, we suggest that the electron-phonon interaction en- 
hances the orbital fluctuations and plays the significant 
role in the realization of the s++-wave superconductivity 
in the realistic parameter region U > U'. In the re- 
cent Raman spectroscopy, it is shown that the electron- 
phonon coupling constant for Aig and Big modes are 
larger {XAig,^Big ^ 0.5)^^ than those predicted by the 
first principle calculations (A ^ 0.21)12,. Then, it is im- 
portant to investigate the effects of the electron-phonon 
interaction on the electronic states and the superconduc- 
tivity based on the microscopic model. 

Recently, Kontani and Onari have investigated the 
5-band Hubbard-Holstein model which includes the 
Coulomb interaction and the electron-phonon interaction 
due to the Big and Eg phonons at the zone center by us- 
ing the RPA and have shown that the Eg phonons drasti- 
cally enhance the orbital fluctuations and the s++-wave 
superconductivity is realized by the orbital fluctuations 
for the realistic values of the electron-phonon coupling^-. 
In ref. 68, they derive the electron-phonon coupling by 
calculating the electrostatic potential variance for Fe-Sd 
electrons from the four surrounding As^~ ions due to the 
oscillations of the Fe atoms assuming that the spatial ex- 
tensions of the Fe-3(i like Wannier functions are small. 
The spatial extensions of the Fe-Sd like Wannier orbitals 
in the 5-band model, however, are very large, e. g., (r'^)- 
(r)2 ~ 5.37A2 for d^2_y2 orbital in LaFeAsO,^"Ii in con- 
trast to their assumption. On the other hand, in the ef- 
fective model which includes both the Fe 3d orbitals and 
the As Ap orbitals, so called d-p model, the spatial ex- 
tensions of the Wannier functions are considered to be 
largely reduced?^"". Therefore, theoretical studies on 
the electron-phonon interaction based on the d-p model, 
are highly desired. 

In the present paper, we investigate the effects of the 



electron-phonon interaction on the electronic states and 
superconductivity based on the two-dimensional 16-band 
d-p model, where the Aig, Big and Eg phonons at the 
zone center are considered. Solving the linearized Eliash- 
berg equation with the pairing interaction obtained by 
using the RPA, we obtain the phase diagram on the pa- 
rameter plane of the Coulomb and electron-phonon in- 
teractions. 

The paper is organized as follows. In Sec. II, we intro- 
duce the the two-dimensional 16-band d-p model coupled 
with the Aig, Big and Eg local phonons and explain the 
formulation of the RPA with the Coulomb and electron- 
phonon interactions in the multi-orbital system. In Sec. 
Ill, we show the numerical results of the charge-orbital 
susceptibility and the gap function for the various values 
of the Coulomb and electron-phonon interactions. The 
linearized Eliashberg equation is solved and we obtain 
the phase diagram. Finally, we summarize the paper in 
Sec IV. 



II. MODEL AND FORMULATION 

Our model Hamiltonian is the two-dimensional 16- 
band d-p model^"— i^i^ coupled with local phonons, 
where 3d orbitals {d3z2_^2, dx2_y2, dxy, dyz, dzx) of two 
Fe atoms (Fe^=^, Fe^=i3) and 4p orbitals {px, Py, Pz) 
of two As atoms are explicitly included. It is noted 
that X, y axes are directed along second nearest Fe-Fe 
bonds. Here, we number the Fe-3(i orbitals as follows: 



^32^ 



^(1), dj.2_y2(2), dxy{3), dyz{4.), dzx{5). 



The total Hamiltonian of the d-p model is given by 



H — Ho + iJii 



H, 



ph 



Hpi^ 



phi 



(1) 



where Hq, i?int, Hph and Hei-ph are the kinetic. 
Coulomb interaction, phonon and electron-phonon inter- 
action parts of the Hamiltonian, respectively. The kinetic 
part of the Hamiltonian is given by the following tight- 
binding Hamiltonian, 

Ho = ^£t4eadiea+^£^plmaPima 

, V^ ,dd >t J.„, I \^ fPP J „. , 

+ H *Z-e.^7n4eaP3ni'^ + f^-c-^ (2) 

where dua is the annihilation operator for Fe-Sd electrons 
with spin a in the orbital i at the site i and Pima is 
the annihilation operator for As-Ap electrons with spin a 
in the orbital m at the site i. In eq. ([2]), the transfer 
integrals tf^^^gg,, tr,j,m,m'> *fj/,m and the atomic energies 
ef, e^ are determined so as to flt both the energy and 
the weights of orbitals for each band obtained from the 
tight-binding approximation to those from the density 
functional calculation for LaFeAsO and are listed in refs. 
64 and 70. The doping concentration x corresponds to 




FIG. 1: (Color online) FSs obtained from the d-p model for 
X = 0.1. The solid and dashed lines show the FSs which have 
mainly dy^jd^x and dx2_y2 orbital character, respectively. 



the number of electrons per unit cell n = 24 + 2a; in the 
present model. The FSs for a; = 0.1 are shown in Fig. 
[T] and we see the two hole FSs (FSl and FS2) and the 
two electron FSs (FSS and FS4) as predicted by the first 
principle calculations^^"^^. 

The Coulomb interaction part of the Hamiltonian is 
given as follows, 

i £ a^a 

+ ^ ^' E E E '^LC' d^u. da. 



2^ E E E dliAi.'d^e'^'d^ea 



(3) 



where U and U' are the intra- and inter-orbital direct 
terms, respectively, and J and J' are the Hund's coupling 
and the pair-transfer, respectively. We assume that the 
relations between Coulomb matrix elements U = U' + 2J 
and J = J' are satisfied throughout the present paper. 

Now we consider the effect of the phonon and the 
electron-phonon interaction parts of the Hamiltonian 
Hph and Hei-ph- By performing the group theoretical 
analysis, it is found that there are 14 kinds of the optical 
phonon modes at the zone center. In the present paper, 
we consider the Aig , Big and Eg phonon modes in which 
As atoms oscillate along the z-axis, Fe atoms oscillate 
along the z-axis and Fe atoms oscillate in the x-y plane, 
respectively (see Fig. [2]). Here and hereafter, we neglect 
the momentum dependence of the electron-phonon cou- 
plings and that of the phonon frequencies for simplicity. 
The resulting phonon and the electron-phonon interac- 
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FIG. 2: (Color online) Schematic figure of the oscillations of 
Fe and As atoms in the (a) Aig, (b) Big and, (c) Eg modes, 
respectively. The small and large spheres denote Fe and As 
atoms, respectively. 



tion parts of the Hamiltonian are given as, 

^P'^ ^^^^sblhs, (4) 

i S 

Hei-ph = E E E E 9fdlAi'M + hsl (5) 

i s 1/1 a 

where his is the annihilation operator for the phonon of 
the mode s {— Aig, Big, Eg and Eg) at the site i, Ws 
is the phonon frequency and gff is the electron-phonon 
coupling. We note that Eg and Eg correspond to the 
oscillation along the x- and y-axis in the Eg mode. As 
following ref. 68, we expand the electrostatic potential 
variance for Fe-Sd electrons from the four surrounding 
As^~ ions due to the oscillations of the Fe atoms for 
the Big and Eg modes in the displacement of the Fe 
atoms up to the first order and expand that in the x, y 
and z coordinates up to the second order. The resulting 
electron-phonon coupling matrix elements of the Big and 
Eg phonons are given as follows. 



Vi9],\^gl]^g%\ 




-V3ghi=gl\=-gf. 


(6) 


9t.=-9%l-V3/2g],l, 


(7) 


9f = 9':\ 


(8) 


g" = (otherwise). 


(9) 



In addition, we also consider the electron-phonon cou- 
pling for the Aig phonon. 



w u X 

9ai„ = 9Aia^e 



(10) 



Within the RPAM^—, the spin susceptibility x^ (?) and 
the charge-orbital susceptibility x'^{q) are given in the 
50 X 50 matrix representation as follows^"— i^, 

riq) = [i^x^°\q)sr'x^°\q), (11) 

x%q) = [i + x^'Hq)C]-'x^°Hq) (12) 

with the noninteracting susceptibility 



iO) a,l3 



Y^ ' ''^ (a) 



T 

'n 



E<".(^)<^.(^ + '^)'(13) 



where a, /3 {—A, B) represent two Fe sites, £ represents 
Fe 3d orbitals, G(fc) = [(ie„ + ^)i - 7?o(fc)]"^ is the 
noninteracting Fe-3d electron Green's function in the 
10x10 matrix representation, fi is the chemical poten- 
tial, Ho{k) is the kinetic part of the Hamiltonian with the 
momentum k in cq. ([2|), k = {k,ien), q ~ ic[,ivm) and 
En = (2n + 1)ttT and v^^ = 2rm:T are the fermionic and 
bosonic Matsubara frequencies, respectively. It is noted 
that when the largest eigenvalue Agpin (Ac-o) of x^^^q)S 
{—x''^Hq)C) reaches unity, the magnetic (charge-orbital) 
instability occurs. 

In the RPA, generally, we need to collect all the ring- 
and ladder-type Feynman diagrams, where the bare ver- 
tices for the spin and charge-orbital susceptibilities S and 
C in eqs. pTt and (IT^ are given by^ 






a, 13 



a,P 



a, 13 



(un^lisi.^ (14) 

25^0Y.9i'''9i'''Ds{iy„r). (15) 



where the Dgi^ifm) = 2uJs/{Vm + "^s) is the local phonon 
Green's function for the mode s and U'^ and C/"^ are the 
bare vertices due to the Coulomb interaction given as 
follows, 

U (U) {a = /3, £i = ^2 = 4 = h) 

U' (-t/' + 2 J) (a = ^, £i = ^3 ^ ^2 = h) 
W iU') = { J {2U' -J) (a = /3, £i = ^2 ^ 4 = U) ■ 
J' (J') {a = /?, ^1 = 4 ^ ^2 = ^3) 

(otherwise) 

(16) 
In eqs. (ITil) and ([TS]), we neglect the ladder terms for 
the phonon-mcdiated interaction. This is valid when the 
condition ujs <C Ep is satisfied, where Ep is the Fermi 
energy, because the ladder terms are proportional to the 
power of uJs/Ep in the weak coupling regime under the 
condition ujs <C Ep^. The vertex corrections includ- 
ing the ladder terms, however, play significant roles in 
the intermediate and strong coupling regime even though 
uJs ^ Ep^. The effect of the ladder terms will be dis- 
cussed in Sec. Ill D. 

The linearized Eliashberg equation is given by 



k' lil2l3lia',P' 



(17) 



where A"^,(A:) is the gap function and V"f ^ ^ [q) is 
the effective pairing interaction for the spin-singlet state. 
Within the RPAH-^^, V^J^i,^^^ (q) is given in the 50 x 50 
matrix, 

V{q) = ^Sr{q)S - lcnq)C + i (5 + C) . (18) 

The linearized Eliashberg equation ([T7|) is solved to ob- 
tain the gap function A"f,{k) with the eigenvalue Asc- 
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FIG. 3: (Color online) Several components of the charge- 
orbital susceptibility x'^{q,0) for U' — and 2g'^/uJo = 0.31 
(a) and those for U' = 1.0 and 2g^ /ujo = 0.45 at J = J' = 0.1 
(b). 



At T — Tc, the largest eigenvalue Age becomes unity. 
In the present paper, we only focus on the case with 
X ~ 0.1, where the superconductivity is observed in 
the 1111 system^. For simplicity, we set T = 0.02eV 
and uJAi„ = ws, — ujp-i = uip2 = loq = 0.02eV 
in the present stud y^^'^^i^^ . As ref. 68, we assume 

and also set 



„25 



„34 
9 pi 



5^2 



„24 _ 44 _ 
9p2 - 5Bi, - 



„35 „ 
9e2 - 

JBlg - 9Alg 



,12 



We use 32 X 32 fc point meshes and 512 Matsubara 
frequencies (— 51l7rT < e„ < SIIttT) in the numerical 
calculations for eqs. pT|) - p^ . and perform the sum- 
mation of the momentum and the frequency in eqs. (jl3[) 
and P?|) by using the fast Fourier transformation (FFT). 
Here and hereafter, we measure the energy in units of cV. 
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FIG. 4: (Color online) The diagonal components of the gap 
function A(fe,i7rr) in the band representation for U' = 0, 
2g^/Loo = 0.31 and J ^ f = 0.1. (a), (b), (c) and (d) cor- 
respond to the 11th, 12th, 13th and 14th bands, respectively. 
The solid and dashed lines represent the FSs and the nodes 
of the gap function, respectively. 
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FIG. 5: (Color onhne) The diagonal components of the gap 
function A(fe,i7rr) in the band representation for U' = 1.0, 
2gV^o = 0.45 and J = J' = 0.1. (a), (b), (c) and (d) cor- 
respond to the 11th, 12th, 13th and 14th bands, respectively. 
The solid and dashed lines represent the FSs and the nodes 
of the gap function, respectively. 
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III. CALCULATED RESULTS 



A. Charge-orbital susceptibility 



As presented in eqs. p4|) and ([T5|) . where we neglect 
the ladder terms, the spin susceptibility is not affected 
by the electron-phonon interaction and is the same as 
the results shown in our previous papers^"— i^. There- 
fore, first, we focus on the effects of the electron-phonon 
interaction on the charge-orbital susceptibility. The ef- 
fects of the ladder terms on the spin and charge-orbital 
susceptibilities will be discussed in Sec. Ill D. 

The several components of the static charge-orbital 
susceptibility x''{q,0) for U' — and 2g^/ujo — 0.31 and 
those for U' = 1.0 and 2g^/ujo = 0.45 at J = J' = 0.1 
are shown in Fig. [3] (a) and (b), respectively, where the 
parameters are chosen to satisfy the condition Age '^ 1- 
It is found that for U' — and 2g-^/uJo = 0.31, the diago- 
nal components oi x'^{q,0), especially [x'^{q,0)]22.22^ ^^^ 
large and have sharp peaks around q ^ (tt, tt) which orig- 
inate from the nesting between the hole FSs and the elec- 
tron FSs (see Fig. [1]), while the off-diagonal components 
are small. When [/' = and 2g^/u}o = 0.31, Ac-o - 0.99 

and the charge susceptibility J2e,e',a,i3{x''{Q:^)]u,e'e' b^" 
comes almost divergent, where the charge fluctuations 
dominate over the orbital fluctuations. We note that the 
charge-fluctuations are enhanced due to the effects of the 
Aig phonon. 

As shown in Fig. |3] (b), for U' = 1.0 and 2g'^/ujQ = 
0.45, the off-diagonal components of x'^{q,0), especially 

[x%q,Q)]24i4 and [x''iQ,0)]tii2 > are large and have 
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FIG. 6: (Color online) The diagonal components of the gap 
function A(fc,i7rT) in the band representation for U' = 1.5, 
2g^/uJo = 0.34 and J = J' = 0.1. (a), (b), (c) and (d) cor- 
respond to the 11th, 12th, 13th and 14th bands, respectively. 
The solid and dashed lines represent the FSs and the nodes 
of the gap function, respectively. 



broad peaks around q = (0,0) and q ^ (tTjTt), while 
the diagonal components are not so large in contrast to 
the case with C/' = and 2g'^/uja = 0.31 (see Fig[3] (b)). 

We note that [x'^IQi 0)]24'.24 and [x'^{q^^)\2iA2 represent 
the transverse orbital fluctuations. The peaks around 
q — (0,0) and q ^ (tt, tt) originate from the scattering 
in the electron FSs and the nesting between the hole FSs 



and the electron FSs, respectively. When U' = 1.0 and 
2g'^/LUQ — 0.45, Ac-o ^ 0.97 and the off-diagonal compo- 
nents [x'^{q)]24 24 ^^'^ [x'^(9)]24 42 ^'"^ largely enhanced, 



^a,|3 



is not 



J24,24 

while the charge susceptibility ^^ ^, ^ gix'^iq)] 
enhanced because of the negative contributions of the or- 
bital susceptibilities such as [x'^('Z,0)]44 55 (see Fig. [3]). 
Then, the orbital fluctuations dominate over the charge 
fluctuations in contrast to the case with U' = and 
2g'^/uJo = 0.31. We note that the off-diagonal com- 
ponents of the orbital susceptibilities [x'^('7i0)]24 24 ^^^ 
[^(9,0)] 

phonon, while [x''(Q', 0)]44^4 and [x''{q,0)]f^% are en- 
hanced due to the effects of the Big phonon. In addition, 
the orbital fluctuations are enhanced also due to the ef- 
fect of the inter orbital Coulomb interaction U' as shown 
in our previous paper—. 



24 42 are enhanced due to the effects of the Eg 



B. Gap function 



By solving the linearized Eliashberg equation, we ob- 
tain the gap function with the lowest matsubara fre- 
quency in the orbital representation A"^ {k,iTrT). Then 
performing the unitary transformation, we obtain the di- 
agonal components of the gap function in the band repre- 
sentation A{k,inT). Figs. Ushows A{k,mT) for f/' = 
and Ig^ jujQ = 0.31. It is found that the pairing symmetry 
is the s+_|_-wave state, where the gap function has no sign 
change on the whole FSs. In this regime, the s++-wave 
superconductivity is mediated by the charge fluctuations 
which is enhanced due to the effects of the A\g phonon 
as mentioned before (see Fig. [3] (a)). 

Figs. [S] shows A{k,iTrT) for U' = 1.0 and 2g'^/ujo = 
0.45. It is found that the pairing symmetry is the s++- 
wave state, where the gap function has no sign change on 
the whole FSs. In this regime, the s+4.-wave supercon- 
ductivity is mediated by the orbital fluctuations which is 
enhanced due to the cooperative effects of the Big, Eg 
phonons and the inter-orbital Coulomb interaction U' as 
mentioned before (see Fig. [3] (b)). 

In contrast to the above two cases, for U' — 1.5 and 
2g'^/u!o — 0.34, the gap function have nodes on the FS4 
as shown in Fig. [Sj Thus, the pairing state is the nodal 
s^-wave state which originate from the coexistence of 
the orbital fluctuations with q = (0, 0) and the spin 
fluctuations with q ^ (tt, tt)^"— i^. When we further 
increase the Coulomb interaction, the spin fluctuations 
dominate over the orbital fluctuations resulting in the 
s^-wave state, where the sign of the gap function be- 
tween the hole FSs and the electron FSs as shown in 
Fig. [ Tp^i^^ . Then, the nodal s±-wave state is observed in 
the crossover region between the s±-wave phase and the 
s+-|--wave phase as shown in Fig. [71 



^ 



^ 



S++ nodal s+ (FS3) nodal s+ (FS4) s+ 



"^ 



_^ ti 



J' 



%-. 
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FIG. 7; (Color online) The schematic figure of the crossover 
between the s++-wave and the nodal s±-wave states and that 
between the nodal s±-wave and the s±-wave states. The solid 
and dashed lines represent the FSs on which the signs of the 
gap functions are positive and negative, respectively. 
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FIG. 8: (Color onfine) (a) The phase diagram on the U'- 
2g^/ujQ plane for J = J' = 0.1 at a; = 0.1, T = 0.02. The 
open and sofid circles represent the s++-wave and s±-wave su- 
perconducting instabilities, respectively. The open triangles 
and squares represent the nodal s±-wave superconducting in- 
stabilities whose nodes are on the FSS and FS4, respectively. 
The dashed, solid, dotted, and dot-dashed Unes show the mag- 
netic order with q ~ (vr, n), orbital order with q ~ (vr, tt), or- 
bital order with q — (0,0), and charge order with q ~ (vr,7r), 
respectively. 



C. Phase diagram 

The phase diagram on U'-2g'^ /loq plane is shown in 
Fig. IHl It is found that the phase diagram includes the 
charge, orbital, and magnetic order and the supercon- 
ductivity, where the charge-orbital, magnetic and super- 
conducting instabilities are determined by the condition 
that Ac-o, Aspin and Asc reach unity as mentioned in Sec. 
II, respectively. 

First, we focus on the ordered phase. As shown in Fig. 
[HI for < [/' < 0.44, the charge order with q ^ (tt, tt) 



takes place at a certain critical value of 2g'^ /ojq. Since 
both the intra- and inter-orbital direct terms U and U' 
suppress the charge fluctuations, the critical values of 
2(7^/ Wo increase with increasing U' . For 0.44 < U' < 1.8, 
the ferro orbital order with q = (0, 0) takes place, while 
the antiferro-like orbital order with q ^ (ti", tt) takes place 
for 1.8 < U' < 1.96. Since the inter-orbital direct term 
U' enhances the orbital fluctuations®^, the critical values 
of 2g'^/uJo decrease with increasing U' . For U' > 1.96, 
the stripe-type AFM order with q ~ (7r,7r) takes place 
as presented in the previous papers^"— i^. 

Now, let us bring our attention to the superconduc- 
tivity. It is found that the pairing symmetry is always 
s-wave and the gap structure sensitively depends on U' 
and 2g'^/uJo. For U' < 1.25, the s++-wave supercon- 
ductivity is realized near the charge and orbital ordered 
phases. The superconductivity near the charge ordered 
phase is mediated by the charge fluctuations character- 
ized by [x'^{q,0)]gg'gg, while that near the orbital ordered 
phase is mediated by the orbital fluctuations character- 
ized by [x''iq,0)]2i%, [x%Q,0)]2iA2^ [x''iQ,0)]ti% and 
[x'^(9jO)]44'55 (see Fig. [3]). On the other hand, for 
U' > 1.6, the s±-wave superconductivity is realized near 
the orbital and magnetic ordered phases, where the spin 
fluctuations are responsible for the superconductivity as 
presented in the previous papers^"— i^. In addition, the 
nodal s±-wave superconductivity is realized between the 
s+^-wave and s±-wave superconducting phases. 

This crossover behavior is naturally explained by the 
relative strength of the orbital fluctuations and the spin 
fluctuations. When the spin fluctuations are not so 
strong, the s++-wave state is realized. With increasing 
U', the spin fluctuations with q ^ (7r,7r) develop, espe- 
cially [x''(q,0)]22' 22- The gap function A^2'^(fe,z7rT) has 
sign change between the hole FSs and the electron FSs, 
where the amplitude of A^2'^(fc,i7rT) around the elec- 
tron FSs is small. On the other hand, Af^{k,iTTT) and 
Af^{k, inT) have no sign change. As a result, the nodes 
appear on the FSS which has mainly orbital 2 character. 
As the spin fluctuations with q '^ (tt, tt) further increase, 
the amplitude of A2-^{k,inT) around the electron FSs 
increase and the nodes appear on the FS4 and finally, for 
U' > 1.6 , the fully-gapped s±-wave state is realized. As 
shown in Fig. [3 the nodes firstly appear on the FSS and 
the position of the nodes smoothly moves to FS4 as U' 
increases. Finally, the nodes on the FSs disappear and 
the s±-wave state is realized. 

We note that the critical value of 2g'^/uJo at which 
the orbital order takes place decrease with increasing 
U' as mentioned above, since U' enhances the orbital 
fluctuations^^. Thus, we stress that the orbital order 
and the orbital fluctuation-mediated superconductivity 
are driven by the cooperative effects of the Coulomb in- 
teraction and the electron-phonon interaction. 



antiadiabatic limit 




U' (eV) 

FIG. 9: (Color online) (a) The phase diagram on the U'- 
2gV^o plane for J = J' = 0.1 at x = 0.1, T = 0.02 in the 
antiadiabatic limit. It is noted that the legends are the same 
as Fig. \8\ 



D. Antiadiabatic limit 

Finally, we discuss the effects of the ladder type dia- 
grams of the electron-phonon interaction which are ne- 
glected in eqs. (fT4|) and ([T5|) . In general, it is difficult 
to include the effects of the ladder terms for the phonon- 
mediated interaction. Therefore, in this subsection, we 
study the effects of the ladder terms by taking the an- 
tiadiabatic limit {ll)s — > oo), where the phonon Green's 
function becomes Dsiivm) — ^ 2/bJs. The resulting bare 
vertices S and C are given by. 



NQ/3 



al3 



(s)l%4.i. = (UX 






28^(,Y^g'f^g'f^lu:,, (19) 

S 

mt,^,,^,^-2b^,Y.^gi^''^^''i^^ 

S 

26^pY.9i'''9i'''/^s. (20) 



Substituting eqs. dH]) and ^ into eqs. ([n]), dH]) and 
([T8| instead of eqs. ^^ and ([15]), we obtain the RPA 
results with including the ladder terms of the electron- 
phonon interaction. Although the antiadiabatic limit is 
not physically relevant, it is still a useful point of the ref- 
erence for getting an overall understanding of the physics 
of the model. We note that, when we apply the antiadi- 
abatic limit to the vertex C given in eq. ([T5|) where the 
ladder terms of the electron-phonon interaction are ne- 
glected, the phase boundaries for the charge, orbital and 
magnetic orders shown in Fig. [5] are unchanged although 
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the superconducting phase boundary is modified due to 
the retardation effect of the electron-phonon interactions 
as mentioned below. 

The phase diagram on U'-2g'^ /ujq plane in the antiadi- 
abatic limit is shown in Fig. [9l As shown in eq. (TTOl) . the 
spin fluctuation enhanced due to the fst term of r.h.s. 
is suppressed by the 2nd term of r.h.s. corresponding to 
the contribution from the ladder terms of the electron- 
phonon interaction. Similarly, the charge-orbital fluctua- 
tion enhanced due to the 2nd term of r.h.s. in eq. ((20|) is 
suppressed by the 3rd term of r.h.s. corresponding to the 
contribution from the ladder terms. Therefore, the mag- 
netic, charge and orbital instabilities are considered to be 
suppressed by the ladder terms of the electron-phonon in- 
teraction. In fact, as shown in Fig. [51 the critical values 
of Ig^ jujQ for the charge and orbital orders with includ- 
ing the ladder terms are larger than those without the 
ladder terms shown in Fig. [5] The critical value of U' 
for the magnetic order with including the ladder terms is 
also larger than that without the ladder terms (see Figs. 
[8]and[9l). 

Remarkably, the s+-|_-wave superconductivity is ob- 
served in a considerably wide parameter region as shown 
in Fig. [S] In the antiadiabatic limit, the bare vertex 
C given in eq. ((20)) is independent of the frequency, 
and then, the effective pairing interaction ^(g) due to 
the charge- orbital fluctuations becomes attractive for a 
wide frequency range. On the other hand, in the case 
with the frequency dependent vertex C given in eq. (|15p , 
the effective pairing interaction V{q) due to the charge- 
orbital fluctuations becomes attractive only for the low 
frequency \vm\ < clIq. As the result, the s+-|_-wave su- 
perconductivity due to the charge-orbital fluctuations is 
observed in a relatively narrow region as shown in Fig. [5] 
for r = Wo = 0.02, where t^(q) is attractive only for Vra 
with TO = as I'l = 27rT > ojq. When the temperature 
is lowered below wq, ^(9) becomes attractive also for Vm 
with |m| > 1, and then, it is expected that the s++- 
wave superconductivity is observed in a wider parameter 
region also for the case with the frequency dependent 
vertex. 



IV. SUMMARY AND DISCUSSION 

In summary, we have investigated the two-dimensional 
16-band d-p model coupled with A\g^ Big and Eg lo- 
cal phonons using the RPA and have obtained the phase 
diagram including the magnetic, charge and orbital or- 
dered phases on the parameter plane of the Coulomb and 
electron-phonon interactions as follows: (1) For weak 
Coulomb interaction, the charge order with q ^ (tt, vr) 
takes place due to the effect of the electron-phonon in- 
teraction with Aig mode. (2) For intermediate Coulomb 
interaction, the orbital order with q ~ (0,0) takes place 
due to the cooperative effects of the Coulomb interac- 
tion and the electron-phonon interaction with Big and 
Eg modes. The orbital order with q ^ (tt, tt) also takes 



place for relatively larger value of the Coulomb interac- 
tion. (3) For strong Coulomb interaction, the stripe-type 
antiferromagnetic order with q ^ (tt, tt) takes place due 
to the effect of the Coulomb interaction. 

Using the effective pairing interaction obtained from 
the RPA, we have also solved the linearized Eliashberg 
equation to obtain the superconducting phase diagram 
including the three types of s-wave pairing as follows: (1) 
Near the charge ordered phase for weak Coulomb inter- 
action, the s+-)_-wave pairing is mediated by the charge 
fluctuations. (2) Near the orbital ordered phase for in- 
termediate Coulomb interaction, the s++-wave pairing is 
mediated by the orbital fluctuations when the spin fluctu- 
ations are not so strong, while the nodal s±-wave pairing 
is mediated by both of the orbital and spin fluctuations 
when the spin fluctuations are rather strong. (3) Near the 
magnetic ordered phase for strong Coulomb interaction, 
the s±-wave pairing is mediated by the spin fluctuations. 

Similar phase diagram including the magnetic and 
orbital ordered phases had recently been obtained by 
Kontani and Onari using the 5-band Hubbard-Holstein 
model®^. They had also discussed the superconductivity 
due to the magnetic and orbital fluctuations and have 
claimed that the s++-wave superconductivity is realized 
near the orbital ordered phase, while the s±-wave super- 
conductivity is realized near the magnetic ordered phase, 
although the detailed superconducting phase diagram 
has not been obtained. In the present study, we have ex- 
plicitly obtained the superconducting phase diagram and 
have found that not only the s++-wave but also the nodal 
s±-wave superconductivity is realized near the orbital or- 
dered phase in contrast to the prediction in ref. 68. In 
addition, the charge order and the charge fluctuation- 
mediated s++-wave superconductivity have been found 
to take place due to the effect of the Aig local phonon 
which was not been considered in ref. 68. In early theo- 
retical studies for the copper oxide superconductors, the 
effect of the Coulomb interaction between the d and p 
electrons Upd was found to enhance the charge fluctu- 
ations which induce the s-wave superconductivit y ^^'^° . 
We have also discussed the effect of Upd on the present 
d-p model for the iron-based superconductors and have 
found that Upd enhances the charge fluctuation-mediated 
s++-wave superconductivity. The detailed results will be 
published in a subsequent paper. 

It seems that both the s±-wave and the s++-wave 
states with full superconducting gaps are consistent with 
various experiments in the iron-based superconductors 
as mentioned in Sec. I, although the sign of the gap 
function has not been directly observed. However, the 
recent theoretical studies of the nonmagnetic impurity 
effects'' - revealed that the Anderson's theorem is violated 
for the s±-wave superconductivity in contrast to the ex- 
perimental results of very weak Tc suppression in Fe site 
substitutio n^^'^^ and neutron irradiation^-. As the impu- 
rity potential due to the Fe-site substitution is considered 
to be diagonal and local in the orbital basis according to 
the first principle calculation®^, it is expected that the 



s+-i--wave state is more robust against the nonmagnetic 
impurity than the s-t-wave state. In the recent ultrasonic 
measurements^!^, a remarkable softening of the elastic 
constant is observed at low temperature down to Tc and 
is well accounted for by Jahn- Teller modes coupled with 
strong orbital fluctuations^. In addition, the weak T- 
dependence of 1/TiT— i2i above Tc in the electron-doped 
compounds is considered to indicate the weak spin fluc- 
tuations. Thus, the s++-wave state due to the orbital 
fluctuations seems to be responsible for the fully gapped 
superconductivity in the iron-based superconductors. 

In BaFc2(Asi_a;P2:)2, the recent field-angle resolved 
specific heaths, and the ARPES measurements^ suggest 
that the superconducting gap function has vertical line 
nodes along the fcz-axis on the electron FSs. This nodal 
superconductivity seems to correspond to the nodal s±- 
wave state obtained in the present study. In the previous 
worka^2i^i22j the similar nodal s±-wave states have been 
obtained in the crossover region between the s±-wave 
phase and the d-wave phase when the different modes 
of the spin fluctuations coexist. This is a striking con- 
trast to the case with the present study where the nodal 
s-l--wave state is realized in the crossover region between 
the s++-wave phase and the s±-wave phase when the 
strong orbital and spin fluctuations coexist. If the fully 
gapped superconductivity widely observed in the iron- 
based superconductors is the s+_|_-wave state, it is natu- 
ral to consider that the nodal superconductivity observed 
in BaFe2(Asi_2:Pa;)2 is the nodal s±-wave state obtained 
in the present study. 

In the present and previous papers^^, we have 
shown that the electron-phonon interaction plays im- 
portant roles for the iron-based superconductors in co- 
operation with the Coulomb interaction. Actually, 
the Raman spectroscopies indicate the large electron- 
phonon interaction^. The large value of the Griineisen 



parameter— and the drastic softening of the elastic 
constan t^^i^^ observed in BaFe2_2;COa;As2 also indicate 
the large electron-lattice coupling. Remarkably, the re- 
cent ultrasonic measurements revealed that the soft- 
ening of the elastic constant C44 is much larger than 
(Cii — Ci2)/2 and continues down to T^^, where the 
temperature dependence of the elastic constant is well 
accounted for by Jahn- Teller modes which couple with 
the orbital fluctuation between dyz and dzx orbitals: 

[x''{l)\lfM ~ [x^ioMi^- Since this type of the orbital 
fluctuation is enhanced due to the electron-phonon inter- 
action with Big mode, we may expect that the effects of 
the Big phonon is most dominant for the elastic soften- 
ing and the superconductivity. In fact, the first principle 



calculation and the Raman spectroscopie: 



,67.78 



suggest 



that the frequency of the Big phonon is lower than the 
Aig and Eg phonons, although the same frequencies are 
assumed in the present study for simplicity. Therefore, 
we need further investigation of the electron-phonon in- 
teraction with including the more realistic effects such 
as the mode dependence of the phonon frequencies and 
the coupling constants and the phonon dispersions which 
have not been considered in this paper. 
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